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Fig. 9. — Grain at the position G is subjected to the flux in the direction given by the 
relative velocity u and its short axis is directed along J. The angle ip is measured in the 
plane II which is perpendicular to the direction of magnetic field H. 
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alignment of prolate grains caused by Alfvenic perturbations, s=-0.5 
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Fig. 8. — The alignment measure oj for prolate grains (g > 0) under Alfvenic perturbations 
(s - 0.5). 
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Fig. 7. — The alignment measure oj for oblate grains (g < 0) under Alfvenic perturbations 
(s = -0.5). 
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Fig. 6. — The alignment measure oj for prolate grains (g > 0) for s > 0. 
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Fig. 5. — 



The alignment measure oj for prolate grains (g > 0) for s < 0. 
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Fig. 4. — The alignment measure oj for flakes (g = 



-0.5) for s > 0. 
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Fig. 3. — The alignment measure oj for oblate grains (g < 0) for s > 0. 
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Fig. 2. — The alignment measure oj for oblate grains (g < 0) for s < 0. 
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Fig. 1. — Zi axis of the external or gas reference frame X1Y1Z1 is directed along the magnetic 
held. The internal or grain frame X 2 Y 2 Z 2 is dehned so as Z 2 coincides with the symmetry 
axis of the spheroid. #1, 6 2} and (p 2 are the polar angles in the above reference frames. 
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Table 1: The analytical expressions for (cos 2 #i) corresponding to different values of grain 
non-sphericity g and external anisotropy s. 





g<0 


5>0 


s < 




V^ arcsin V~T+7 1 


s arc tan-* /., , s ^, — s 


s arctanh A /— 7-^; — s 


s > 


\f— oarcsinh. / 77- 


v^ arcsinh A/rfi' 1 


s arctanh A /— . , s ^, — s 


s arctan-v /-r-r-^-. — s 

V l+s+9 
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-oarcsinh, /-T— i 
(cos ^ ^ = _^ _^+£_ _ i ? (B15) 

5 arctanh, /— 5 



V^arcsinh^ x 

(cos ^) = == - -, (B16) 

5 arctan, ' — — 



for 5 > 0, and g < 



for 5 > 0, and g > 



l+s+g 

which cover all the cases. 



C. Alignment of rotationally hot grains 

We have seen in the main body of the paper that if diffusion J dominates its "leaps", 
the density of J in the space of angular coordinates for an ensemble of grains is proportional 
to S}J 2 . Therefore 

fM) = C~ 1 (C17) 

Jo cos^a 

where 

c= r di ,r^, ,ci8) 

Jo Jo cos 1 '^ 

provides an opportunity to obtain the Rayleigh reduction factor for any <f> through numerical 
integration in Eq(C18) 

However in the most important case corresponding to the alignment under Alfvenic 
perturbations, the calculations can be done analytically. Indeed, cos a = sin <p cos ?/> and 
according to Gradshtein & Ryzhik (1965, 3.621(f)), 

r-vr/2 







sin 3 / 2 " 1 xdx = 2 3/2 " 2 5(3/4, 3/4), (CI9) 



where B(x } y) is beta function. Similarly, 

sin 3 / 2 " 1 x cos 2 xdx = 2 3/2 " 2 5(3/4, 3/4) = — , (C20) 

Jo ' 3^25(9/4,1/4) V ; 

where at first cos 2 x was expressed through cos 2x and then Gradshtein & Ryzhik (1965, 
3.631(8)) was consulted again. Expressing Beta functions through Gamma functions: 
5(3/4,3/4) = T 2 (3/4)/r(3/2) and 5(9/4,1/4) = r(I/4)r(9/4)/r(5/2) and substituting 
the corresponding values of T it is possible to obtain 

(7^0.04, (C2I) 

which is much less than the alignment when "leaps" dominate. 
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where the above identity is taken into account. The integral can be calculated [see 
Gradshtein & Ryzhik 1965, 2.271(4)] to give 

dx 



° +g) + SX 2 

l ln ^±^±i (a>0) 



+g 

(B8) 



i 



■arcsm 



where for uniformity one can also use inverse hyperbolic function for s < 0, namely 



1 / s 
u = — =arcsinh,/ . (B9) 



To obtain C(s,g), one has to calculate 

which can be done by substituting u = x 2 + s. The corresponding integral can be easily 
evaluated [Gradshtein & Ryzhik 1965, 2.224(5)] to give 

12={ (Bll) 
|=arctany3^ (sg > 0). 



The function C(s,g) is equal to i 2 1 . Using the inverse hyperbolic function the expression 
for sg < is equal 

1 / —sg 
to = — arctanhw , (B12 

which enables one to find C(s,g). 

As a result, one obtains for s < 0, and g < 



-flarcsm./-— i 

(cos 2 ^) = - Y=i+£._i (B13) 

s arctan,/^^ s 



for 5 < 0, and g > 



./Tfarcsin, /— -7— 1 

(cos^ ^) = - i, (B14) 

5 arctanh, /— ■§ 
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b 00 = c(l + s(l-xi)+g(l-x 2 )+gs(l-2xi)(l-x 2 )) 

hi = -^(1 + 5^ + 1(1 + ^)^(1-2^)) 

b 22 = j^(l+gx 2 2 + ^(l + x 2 1 )gsx 2 2 (l-2x 2 1 )) 

boi = -7T - xl)s(l + g(l - x 2 2 )) 

bo2 = x 2 (l - x 2 )g(l + 5(1 - x\)) 

W2 = —xix 2 sg(l - xi)(l - x 2 )), 



where 



c = l«P 2 ) - (Pl))b 2 

* = -\((p 2 ) - Hp 2 M(p 2 ) - (pI))- 1 

a 2 -b 2 

9 = 

and p and p z are the values of the momentum and its Zi projection transferred in a collision. 



B. Computation of integrals 



This double integral and the one in Eq. (14) can be integrated over 6 2 to give 

(cos 2 0^ = (s,g) f / dX (B4) 
Jo (1 + sx z )y/L + $x z + g 

and 

c "-"jr ( i +M »M +M » +g Bi - (B5 » 

where an evident substitution x = cos 61 is used. Writing 

A(l + sx 2 ) + B = x 2 (B6) 

it is possible to obtain A = 1/s and B = —1/s. Therefore 

1 /"l dr 1 

(cos 2 1 ) = C(s,g)- --, (B7) 

s Jo y 1 + sx z -\- g s 
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A. Coefficients of the Fokker-Planck equation 



The coefficients a 4 - and bik of the Fokker-Planck equation are determined by the change 
of the angular momentum of the grain due to its collision with an atom: 

/ 1 dAx x 1 A dAx x A \ /A , 

a = ( - A x 2 — — - + - A <^i— ^ — - - Ax 1 ) Al 
\2 ox 2 2 ocpi I 

1 1 9 A x m 1 d Ax m 

ZA X\ 1 ZA X r 



\2 dxi 2 9x m 

1 A d A x m A \ , 4 . 

+ -Ay m -Ax m (A2) 

Z O ZA (,C m / 

= (Ai, A , (A3) 

where m = 1, 2, z = 0, 1, 2, and the angular brackets denote averaging over the impacts 
of atoms over the grain surface. The quantities Ai, and A(fi and the corresponding a 4 - and 
bik were calculated in Dolginov & Mytrophanov [1976] using the following equations: 

Axj = A(ej ■ J) - (ej • J) A J 
Ax = j-AJ 

A^- = (( ej -xj)-AJ(J(l-a?))-\ 

where j = 1,2, ei = is a unit vector along the magnetic field, e 2 = j|j is a unit vector 
along Zi-axis of the grain, j = ^ is a unit vector along J, and AJ = r X p. 



To find coefficients (k = 0, 1, 2) and bik (h k = 0, 1, 2) one needs to substitute A 



J5 



Ax , Aifj into Eqs. (Al), (A2), and (A3) and perform the necessary averaging. Apart from 
averaging over the part of surface exposed to the flux, one has to average over the angles of 
precession of J around e 2 and of J around m. As a result, one gets for an axially symmetric 
ellipsoidal grain [Dolginov & Mytrophanov 1976]: 

a = j(2 + s(l + x 2 1 )+g(l + xl) + 2gs(x 2 1 -x 2 2 )) 

ai = -^(-2 + s(l -3x1) -g(l + x 2 2 ) 

+ 2gs{l -2x1-2x1 + 3x1x1)) 

a 2 = — x\{-2 - s(l - x\) + g(l + 3x1) 

+ 2gs{l -2x1-2x1 + 3x1x1)) 
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with their axes of major inertia perpendicular to the direction towards the UV source. For 
th 3> td } t x which is ~ tL becomes again proportional to S~* and a can be computed in 
accordance with Eq.(41). The mechanism will tend to align long axis of grains parallel as 
well as perpendicular to magnetic field lines depending on the angle between the direction 
of magnetic field and that towards the UV source. For typical ISM parameters, one needs 
to account for paramagnetic relaxation if 3> td and therefore the resulting alignment 

(2) 

can be estimated according to Eq.(44), if one denotes <Tj the measure of paramagnetic 
alignment. However, a detailed discussion of this radiation driven mechanism is far beyond 
the scope of our present paper. 

Everywhere above we disregarded paramagnetic relaxation of the suprathermally 
rotating grains. This is justifiable if we adopt the "standard values" of the ISM magnetic 
field and gaseous density and assume < td (see Spitzer 1978). If it is not the case, 
one may use Eq. (44) to estimate the alignment for the joint action of mechanical and 
paramagnetic processes. 

7. Conclusions 

We have shown that suprathermal rotation caused by H 2 formation does not prevent 
grains from being aligned mechanically. The alignment arizes from both gas depositing 
momentum with grain during crossovers and due to the change of the mean time back to 
crossover. These two processes act in the same direction and tend to minimize the cross 
section of grain interaction with a gaseous flux; the second insures that grains exhibit 
alignment even if the random torques are dominated by H 2 formation. Alternatively, if the 
gaseous bombardment dominates random torques, the first process dominates. Our results 
show that the alignment is efficient for oblate grains subjected to Alfvenic perturbations or 
radiative fluxes. They also indicate that the alignment caused by Alfvenic waves can be 
widely spread. 

This work was initiated by Bruce Draine's comments, but would not be possible if 
not for encouragement by Ethan Vishniac. I am grateful to Russell Kulsrud for valuable 
comments on grain charge, Alyssa Goodman and Phil Myers for illuminating discussions on 
observational data, and to David Williams for explaining me the most subtle issues of grain 
chemistry. The manuscript was revised after my visit to Princeton University Observatory 
where I got much from stimulating discussions, especially, with Bruce Draine and Lyman 
Spitzer. It is a pleasant debt to thank Jeremiah Ostriker for arranging financial support 
for this visit. This paper also owes much to my communications with John Mathis, Martin 
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that PAH might be aligned by the ambipolar diffusion and this may produce polarization 
in the associated emission features. 

Everywhere above crossover and spin-up alignments were discussed separately. In fact, 
they act together and we may estimate the resulting Rayleigh reduction factor. In general, 
if the component of J along magnetic field is J Z} the alignment due to a mechanism acting 
alone results in the ratio X{ = , T x 2 ( '\ = , T y 2 ( '\ and the measure of alignment for this particular 

\ J z(i)l \ J z(i)l 

process can be estimated as follows 

^"-f^—--). (42) 
3 2 \2xi + 13/ v ; 

If m independent processes act simultaneously, the corresponding ratio 

determines c^K This gives a way of expressing through 

For example, if the crossover alignment corresponds to and the alignment in the 
sequence of crossovers corresponds to <Tj their joint action corresponds to 

J ~ 1 + 2^*? ■ ' ' 

Thus, the measure of alignment for flakes subjected to Alfvenic perturbations can be 
estimated 0.61 if the two mechanisms act together. 

If, one of the processes dominates and, for instance, a^P ^> <Tj 2 ^ it is easy to see that 

-i D -4 1) + 4 2) (l-4 1) )- (45) 



If recoils from the H 2 molecules formed over grain surfaces dominate the random 
torques, the alignment in the course of crossovers is suppressed and grains can be 
mechanically aligned only by the second mechanism. This is likely to be true for relative 
velocities u just above the sonic ones or if the accommodation coefficient is greater than 
we assumed. Alternatively, if > tj, and, due to some reason, is independent of the 
accretion rate, it is the first mechanism that operates. 

In fact, our study also indicates the existence of another type of alignment mechanism. 
Namely, if is controlled by photodesorbtion and ^ td, the distribution of J will 
become anisotropic even in the absence of any mechanical flows; grains will tend to align 
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dichroic absorption, the alignment is dominated by the difference in cross section, that is 
studied above. As soon as becomes comparable to the Purcell & Spitzer effect becomes 
negligible even for spherical grains. 

However, the fact, that in our simplified model we may ignore "localized" adsorption - 
desorption events does not mean they do not deserve a further study. Note, that for most 
favorable conditions, e.g. when the flow is parallel to magnetic field lines, such "localized" 
events may provide alignment measure of J for spherical grains of the order of 5/17 ~ 1/3. 
This estimate is greater that that in Purcell & Spitzer (1971) as we considered quasi-regular 
torques as opposed to random torques in the latter study. However, a more detailed 
discussion of the consequences of "localized" adsorption - desorption events is beyond the 
scope of the present paper. 

Above grains were approximated by discs. Another extreme corresponds to needles. 
It is easy to show that the measure of alignment for them tends to zero for the fluxes 
perpendicular to magnetic field lines and approaches maximum of the order of 0.1 for fluxes 
along magnetic field lines. This apparent inefficiency of alignment for prolate grains as 
compared with oblate ones stems from a marginal difference in the time averaged grain - 
gas cross sections for different orientations of rotating prolate grain. Therefore in major 
cases it is possible to disregard the contribution from mechanically aligned suprathermal 
prolate grains. Note that there is an observational evidence that the aligned grains are 
typically oblate (see Hindelbrand 1988). 

6. Discussion 

We have seen that suprathermal grains can be aligned by a supersonic gaseous flux. It 
is also shown, that both mechanisms discussed tend to minimize the cross section of the 
flux interaction with the suprathermally rotating non-spherical grain. 11 

As the alignment depends on the relative gas-grain velocity (see section 3), one may 
predict that larger grains in diffuse clouds should be more efficiently aligned by Alfvenic 
waves, as compared to smaller ones. This is a tendency that corresponds to observations 
(Mathis 1979). However, drift velocities caused by ambipolar diffusion do not decrease for 
small grains. Therefore if ambipolar diffusion is responsible for alignment of grains in dense 
clouds, small grains should be also well aligned. If small grains are non-spherical, this trend 
may be detectable. For instance, it was suggested by B. Draine (private communication) 



11 Note, that the Gold mechanism acts in the same way. 
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-i 



2^ dtp 



(39) 



cos a 



where 



/'27Y /' 

C = dip 
Jo Jo 



>7r sin (fid(fi 



(40) 



cos a 



provides an opportunity to obtain the Rayleigh reduction factor for any <f> through numerical 
integration: 



However, for the most important case, corresponding to Alfvenic perturbations, all 
the calculations can be done analytically and a = 0.25. Such an alignment, although less 
efficient than the one we discussed in details in section 3 cannot be ignored. For streaming 
along magnetic field lines a = — 0.5. A joint action of the two mechanisms is discussed 
below. 

Note, that the alignment caused by Alfvenic perturbations is marginal for grains 
rotating due to cosmic ray bombardment (i.e. a = 0.04). This does not mean, nevertheless, 
that the mechanism is not efficient at all. For instance, if the flow is directed along magnetic 
lines, it is possible to show that the corresponding measure of alignment a is equal to 0.2. 
However, it is beyond the scope of the present paper to discuss this effect. 

It worth noting, that the possibility of alignment of grains subjected both to supersonic 
drift and to cosmic ray bombardment was mentioned back in Salpeter & Wickramasinghe 
(1969). However, the authors believed that the alignment would be orthogonal to that 
attainable through the Gold mechanism. On the contrary, we have proved that for 
non-spherical grains the two mechanisms act in the same direction. This was not found in 
Purcell & Spitzer (1971) as there spherical grains were studied, while the effect we speak 
here is caused by grain non-sphericity. The residual marginal alignment obtained in the 
latter paper arizes from the assumed peculiar interaction of atoms with a grain. Namely, in 
the model adopted in Purcell & Spitzer (1971) colliding atoms were not allowed to diffuse 
over grain surface. This poses an interesting question to what extend our results depend on 
the assumption of atoms being adsorbed by the surface on the collision. For spherical grains 
Purcell & Spitzer found that the gaseous friction for the rotational axis perpendicular to 
the flow exceeds 1.5 times that for the rotational axis parallel to the flow. This friction 
influences tj, and therefore for <C td the effect of atoms evaporating from the same spots 
over that they hit the grain surface may act in the opposite direction as compared to the 
effects of crossover difference studied above. However, it is easy to see that for the Alfvenic 
perturbations the Purcell & Spitzer effect is suppressed due to grain precession. Moreover, 
for sufficiently non-spherical grains, which are most interesting as the media responsible for 




(41) 
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In our two dimensional space of angular coordinates, the time intervals At 3 on average 
are proportional to t x . The damping time tj that enters the expression of t x is inversely 
proportional to the rate at which atoms arrive to the grain surface (see Eq.(2)); it also 
natural to assume that t^ varies in the same proportion. Indeed, the time t^ is proportional 
to the time of colliding with N\ heavy atoms which will poison at least half of active sites; 
the rate of heavy atoms arriving to the surface is proportional to S n . Therefore tji, and 
tj x , which are, respectively, the life-time of Purcell's rockets and the mean time between 
crossovers for the position j, are inversely proportional to the corresponding cross section 
Sj n ; the coefficient of proportionality is the problem of normalization. Thus a statistical 
description of the alignment of a suprathermally rotating grain subjected to a supersonic 
flow becomes a matter of determining how grain-gas cross section changes in the course of 
grain precession in the ambient magnetic field. This problem can be solved for grains of 
arbitrary shapes. 

For diffusion of J when t e \ < t e \v , it is possible to see that J scales as ^Jtl ~ S~ x l 2 (see 
Eq (30)), while the deviation AJ scales as a square root of the number of torque events, 
i.e. ~ S^ 2 . Therefore the time At; and tj scales again as S' 1 . This universality of scaling 
for "leaps" and "diffusion" regimes is a consequence of the fact, that both A J and J are 
controlled by the intensity of the flux. If we assume that A J is controlled by some other 
process, e.g. by cosmic ray bombardment (see Salpeter & Wickramasinghe 1969, Sorrell 
1995), while J is limited by friction caused by a supersonic flux, tj would scale as S' 1 ^ 2 . 
We briefly discuss this alignment in Appendix C. 

To simplify our treatment, while elucidating the nature of the effects, consider grains 
approximated by infinitely thin discs; a quantitative treatment of alignment for other 
shapes will be given elsewhere. An advantage of using a disc, rather than other shape is 
that a straightforward expression is available for S n : 



where a is the angle between the direction of the flux defined as the direction of u and that 
of J (J is perpendicular to the disc plane as a result of the Barnett relaxation). If J and u 
make, respectively, angles <p and <f> with H and the angle between the planes uH and JH is 
if) it is obvious from the spherical trigonometry (see Fig 9) that 



The distribution function / should be averaged over ip as this is the azimuthal angle of 
J that changes in the course of precession about H direction. Therefore the distribution 
function 



S n = irr 2 cos a, 



(37) 



cos a = cos <f> cos <p + sin <f> sin <p cos ip. 
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For an ordinary grain td/t e i is of the order 10 9 , which is of the order of the atomic to grain 
mass ratio, while v does not exceed 10 6 even for the surfaces completely covered by active 
sites. Therefore one may expand the square root in Eq (34) and to obtain, that for the 
diffusion to dominate J dynamics either 

t L £ U-^- (35) 

Vt e l 

or 

t L £ Uv. (36) 

The condition given by Eq (35) is not likely to be satisfied for the typical ISM 
conditions (see Spitzer & McGlynn 1979). Moreover this situation is not interesting from 
the point of view of mechanical alignment, as the expected paramagnetic alignment is 
nearly complete for such long t^. 

The condition given by Eq (36) is, in fact, a criterion for the existence of quasi-regular 
torques 10 . Any H atom arriving at the grain surface has nearly equal chances to form 
an H 2 molecule over any of v active sites. As these sites are expected to be distributed 
randomly over grain surface, the grain experiences random torque over the interval less 
than t P jv. Such short are expected for grains less than a critical size (Lazarian 1995a, b). 
For grains smaller than this size, every oxygen atom has good chances to poison an active 
site. As the poisoning time of v active sites is inversely proportional to the oxygen flux, the 
above criterion (see Eq (36)) for / < l cr is equivalent nfj^i/n < 1, where n and njj are 
the concentrations of atomic oxygen and hydrogen respectively. Therefore, it is obvious, 
that for reasonable values of the accommodation coefficient 71, e.g. 71 = 0.2, diffusion 
can dominate leaps only for the core regions of molecular clouds. Above it was implicitly 
assumed that in spite of rapid poisoning there exist an efficient mechanism of creating new 
sites. In fact, H 2 formation is likely to be suppressed for I < l cr . In short, for the majority 
of cases disorientation during crossovers dominates the dynamics of J. 

To find the angular distribution of J, we need to consider a long it — > 00) sequence of 
individual "leaps". We assume that in the course of a sequence of "leaps" J has an equal 
probability of obtaining any direction within [0,27r]. In this case, if J spends time tj in a 
particular volume of the phase space, which is a sum of the At; intervals that it spends any 
time on entering the volume, lim^oo ^ gives the time averaged probability of J entering 
the volume. According to the ergodic hypothesis this average coincides with an ensemble 
one. This way of reasoning is easy to generalize to a continues distribution. 



10 The diffusion can be prevalent over leaps in spite of the existence of quasi-regular torques, 
if the random torques are dominated by gaseous bombardment. 
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becomes not important as compared with the overall angular momentum of the grain, it is 
the rate at which atoms arrive to the grain surface that matters. This rate influences both 
the diffusion of the J and the frequency of crossovers; we will show further on that either of 
these effects is capable to produce an anisotropic distribution of grain axes. 

Due to stochastic torques the direction of J is subjected to diffusion, the characteristic 
time of which can be estimated as the time during which the mean square deviation a/AJ| 



becomes of the order of y (J 2 ). This diffusion is a random walk process and therefore: 

AJ| « AJ e 2 ^, (29) 

where A J P j is the angular momentum deposited in an individual elementary torque event 
and t e i is the mean time between the torque events. As, according to Eq (3), the mean 
squared angular momentum is 



LP) « v-V 2 AJ el - J—^— (30) 
the characteristic time of diffusion is 

tel td + *L 

In general, both AJ e ; and t e \ in Eqs. (29) and (30) are different. For instance, if 
stochastic torques caused by atomic bombardment dominate, AJ e ; and t~l in Eq (29) 
are, respectively, angular momentum and frequency associated with atomic impacts; those 
quantities can differ several times from ones entering Eq (30). However, for our order of 
magnitude estimates this should not be very important. 

Consider now crossovers. From the point of view of J dynamics in angular coordinates 
they are equivalent to leaps; the corresponding time of disorientation as a result of such 
leaps is 

tleap ~ 7 7 T^TT , (32) 

arccos(exp(— r )) 

where we remind the reader, that F is the disorientation parameter. For sufficiently 
small g 7r/2, which corresponds to a complete disorientation during a 

crossover. The time t x was estimated in Purcell (1979) as 

t x = 1.3(t L + t d ). (33) 
It is easy to see that tdijf equals to ti eap if 



t L = -t d + d — ± d — J - + 1. 34 
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ambipolar diffusion in the context of mechanical alignment of grains was first mentioned in 
Whittet (1992) with the reference to a forthcoming paper by Roberge & Hanany. 

According to Draine, Roberge & Dalgarno (1983), grain velocity in respect to neutral 
gas is 



where v n and V{ are the velocities of neutral and ionized components respectively, while u c 
is the grain gyrofrequency and t m is the time that takes a grain to collide with gaseous 
atoms of the net mass equal to that of the grain. According to fig. 1 in Pilipp et al. (1990), 
grain-neutral velocities may be highly supersonic for sufficiently strong shocks. 9 

Hereafter while discussing alignment under Alfvenic perturbations we will bear in mind 
both supersonic grain drifts in highly ionized media due to Alfvenic waves as well as grain 
drift due to magnetohydrodynamic shocks in weakly ionized media. 

The measure of alignment for oblate and prolate grains under Alfvenic perturbations 
is shown in Fig 7 and Fig. 8 respectively. It is obvious that for prolate grains the alignment 
is negligible. At the same time a comparison between fig. 5a in Lazarian (1994a) and Fig. 7 
indicates that although the alignment for flakes is less efficient for suprathermally rotating 
grains as compared with thermally rotating grains, for grains of moderate oblateness, both 
processes deliver a comparable degree of alignment, which is of the order of 20%. 



Suprathermal grains are not sensitive to the angular momentum deposited by chaotic 
gaseous bombardment during spin-ups. However, this does not mean that gaseous flux does 
not cause alignment of suprathermal grains apart from relatively short crossover intervals. 

We will show below that the gaseous flux influences the distribution of the angular 
momentum for an ensemble of non-spherical suprathermal grains. The condition of 
non-sphericity is essential for the mechanism under study, as we will show that it is the 
difference in cross section of the flux - grain interaction that produces the anisotropy. For 
spherical grains, which are the favorite object for theoretical studies, the effect vanishes. 
In other words, as the angular momentum deposited with the grain by the gaseous flux 



9 Note, that if the charge density of grains is a significant fraction of that of the ions and 
electrons the grain drag is independent of the parameter u c t m (Nakano & Umebayashi 1980). 




(28) 



5. Alignment in the sequence of crossovers 
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minimize their cross-section. Therefore a further decrease of the upper limit of is 
expected for non-spherical grains. 

In short, for grains with negligible charge, supersonic drift must be widely spread in the 
ISM for grains of radii above 10~ 6 cm. Charging of grains modifies this conclusion. In fact, 
the charge issue is a subtle topic as the two main causes, i.e. collisions and photoelectric 
emission, result in charges of opposite signs and the grain charge in diffuse clouds depends 
on the interplay of these two processes. To start with, consider the charge caused by a 
disparity in the electron and ion collisional rates. According to Spitzer (1978) the grain 
potential energy eU becomes of the order of — 2.5&T, where T is the ion kinetic temperature. 
A grain with such a charge moves about magnetic field lines with a cyclotronic frequency 

UaB 2.5kTaB , , 

u c = . (27) 

m g c em g c 

For p g = 3 g cm" 3 , a = 10" 5 cm, T = 80 K and B = 3 x 10" 6 G Eq. (27) provides 
lo c 4.6 X 10~ 12 s -1 , which is of the same order that u max . The ratio p 2 = ^c/^A is another 
important parameter of the theory for // 2 < 1 the drift velocity scales approximately as 
v p 2 - This limits the sizes of grains that move supersonically. In fact, only grains with 
radii above 10~ 5 cm are likely to drift supersonically as a result Alfvenic oscillations, if the 
values lo c and u max are given as above. Incidentally this corresponds to the size distribution 
of grains that cause polarization (see Kim & Martin 1994, 1995). If the charge is reduced 
by photoelectric emission or u max is greater, the supersonic drift can persist in the diffuse 
clouds for smaller grains. 

In short, in order to provide supersonic drift of grains of the size ~ 10~ 5 cm 8 , high 
frequency Alfvenic waves should have large amplitudes. This does not contradict to our 
present day knowledge of Alfvenic waves, which "evolve by steepening until ion-neutral 
collisions damp them" (McKee et al 1993) and cannot be dismissed by considering the 
energy balance within galactic interstellar matter either. However more elaborate study is 
needed before any definite far-reaching conclusion can be made. 

Grain charge provides a possibility for alignment through supersonic ambipolar 
diffusion. This process cannot be ubiquitous diffuse clouds as an assumption of supersonic 
ambipolar diffusion over large scales entails energy dissipation well in excess of what 
supernova explosions can inject into the ISM turbulence. However, this may be an option 
for particular regions subjected, for instance, to MHD shocks. Note, that the importance of 



8 If small grains drift with supersonic velocities, then a fortiori larger grains should drift 
supersonically. 
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gas-grain velocity, which is the difference between the atomic v a and grain v g velocities, i.e. 
u = v g — v a . For harmonic perturbations v a = ^osincu^t, and the equation ol motion can be 
written as 

T" + ViP\ -f + P 2 = ~ cos r, (24) 
dr v Vq 

where r = tu>A, and p = u/v 0} while p\ = m ™ S ^"° ~ ^p;^" 6 ? where g g is the grain density, 
a is the grain radius, and loa is the Alfvenic frequency and we remind the reader, that m g 
and m are, respectively, the masses of a grain and an atom. For subsonic motions, the 
solution is 

u= , = sm(io A t + <f) sh ), (25) 



>l + p 2 

where tan (f> s h = p^. This solution reflects the most essential features of the generalized 
problem, particularly, a decrease of the amplitude of u with fj^. The lower limit for this 
coefficient can be expressed as (Lazarian 1994a) 

/^r-c 1 ^^, (26) 

where is the ionization ratio, m a is atomic mass, {(t^vj) 1.5 • 10~ 9 cm 3 s -1 is the 
collision rate coefficient (Nakano 1984). To obtain the above estimate we used the maximum 
frequency of Alfvenic waves that still move ionized and neutral components together 
2i^n(aTVj) (McKee et al 1993). 

Substituting = 10~ 4 , n = 15 cm -3 , a = 10~ 5 cm and p g = 3 g cm 37 , one obtains 
fi™ m 2 • 10~ 2 <C 1, which indicates that supersonic Alfvenic waves should produce 
supersonic grain drift. The corresponding cut-off frequency 4.5 x 10" 12 s~\ 

Photoionization can increase the ionization ratio (McKee 1989) and data in Myers 
& Khersonsky (1995) indicate, that the ionization ratio is high up to the densities 
corresponding to dark clouds. 

For our estimates above we have assumed that a grain is a sphere. Our computations 
above show that under gaseous bombardment grains tend to align in such a way as to 



velocities are subsonic, one should use the total surface of grains, but the flux becomes 
0.25gv s . For a sphere, the force is well approximated by ira 2 gu^Ju 2 + v 2 s for both supersonic 
and subsonic motions. 

6 Note, that ji\ is proportional to the ratio of the Alfvenic wave period to the the time tj. 

7 Following Spitzer (1978) we assume that p g = 3 g cm 3 for a = 10~ 5 cm and decreases to 
1 g cm 3 for grains of radius 5 X 10~ 5 cm. 
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4. 



Alfvenic perturbations &c supersonic grain motions 



Although above we have discussed alignment for various directions of the grain 
drift velocity, not all of the corresponding parameters s appear to be equally important 
when we study grain alignment within the ISM. One of the "chosen" parameters s — > oo 
corresponds to the motion of charged grains under the radiation pressure, when this motion 
is constrained by magnetic field. Estimates in Purcell (1969) show that such a supersonic 
motion is important in the vicinity of stars, but cannot persist over an appreciable part of 
the ISM. Therefore a drift under the influence of the perturbed magnetic field will be our 
major concern below. 

First, we study if grains in diffuse clouds decouple from the moving ionized gas at lower 
frequencies than ions decouple from neutrals and whether this results in the supersonic drift 
for any conceivable values of the ISM parameters. 

It is generally accepted that the Alfvenic supersonic turbulence dominate the dynamics 
of random motions within the ISM (Aron & Max 1975, Myers 1985, Falgarone & Puget 
1986, Elmegreen 1990, Heiles et al. 1992) and waves of high frequency naturally arise 
due to non-linear steepening of linear polarized and unpolarized waves as they propagate 
(Elmegreen 1992, McKee et al. 1993). Energy dissipation that accompanies these supersonic 
motions is believed to be one of the major heating mechanisms for the intercloud (Ferriere 
et al. 1988) and for interclump (McKee 1989) gas. Here we want to know whether these 
motions can cause alignment. 

If the Alfvenic velocity exceeds that of sound, it is possible to show that the motions 
of ionized gas are localized mainly in the plane perpendicular to the magnetic field lines 
(see Alfven & Falthmmar 1963) and this corresponds to s ~ —1/2. Further on, for 
the sake of simplicity, we will assume that for Alfvenic perturbations the drift velocity 
u is perpendicular to H. This is a good approximation for standing Alfvenic waves. 
For travelling Alfvenic waves, in general one needs to perform averaging over angles 
[7r/2,7r/2 — arctan B±/B z ]. However we will disregard this within our simplified approach. 

Both a mechanical force due to gaseous bombardment and an electromagnetic one act 
upon a charged grain subjected to the Alfvenic waves. We discuss at first a mechanical 
force. This force applied to a grain can be expressed as (Kwok 1975) 



where Sd ~ 7ra 2 is the grain cross-section 5 , v s is the sound speed, and u is the relative 




(23) 



5 This cross section should be used when the relative velocities are supersonic. If the 
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of angular momenta and according to Eq. (15) oj = a = 0. The absence of alignment for 
s = is evident from all the figures presented in the paper. As s tends to infinity, (cos 2 0\) 
tends to zero and therefore oj = o tends to —0.5 (see also Fig 4). High values of s that are 
necessary for the alignment are easily attainable when charged grains stream along field 
lines. 

Proceeding with our discussion we consider right upper corner of Table 1 which depicts 
alignment of prolate grains tending to drift at right angles towards magnetic field lines 
(i.e. s < and g > 0). The corresponding oj is drawn in Fig. 5. It is easy to find that 
for needles, i.e. g — > oo, (cos 2 #i) 0.38, which means that oj 0.08. Note, that the 
corresponding Rayleigh reduction factor a —0.04. The alignment measure increases for 
less prolate grains and reaches 



for spherical grains. In short, alignment of prolate grains for s < is not efficient. 

The right lower corner of Table 1 presents alignment of prolate grains tending to drift 
along magnetic field lines (i.e. g > 0, s > 0). It is apparent from Fig. 6 that the alignment 
is efficient for large s. The measure for J vectors tends to —0.5 which corresponds to a value 
of 0.25 in terms of the measure of axis alignment. Therefore streaming along magnetic field 
lines should be efficient in aligning prolate grains with their axes along the lines. 

To summarize, the alignment discussed above is similar to the Gold one in the sense 
that it is caused by the momentum deposited by a gaseous flux with the grain. The 
difference is that the alignment proceeds in short periods of crossovers, when grain is 
susceptible to the gaseous bombardment. This entails the two stage alignment process 
discussed above. 

All the way above it was assumed that the velocity of gas-grain drift substantially 
exceed the thermal one. However, it is possible to estimate the influence of thermal motions 
of atoms. Indeed, in the grain reference frame, these motions are seen as deviations of 
individual atoms from the mean flux direction. The angle of this deviation can be estimated 

as 



and therefore it is possible to obtain the measure of alignment through averaging our results 
over <j. In terms of s parameter used above, the averaging over <j for ? < 1 corresponds to 
averaging over 



which gives an opportunity to use the graphs above to account for a finite gas - grain drift. 




(20) 



q arcsm 







- 10 - 



of J in the grain frame of reference. Therefore according to Eq. (16) we assume a uniform 
distribution of J the internal frame of reference as compared to using 8 function distribution 
for 6 2 assumed in Lazarian (1994a). The calculations relevant to obtaining the analytical 
expressions for (cos 2 0\) are performed in Appendix B and the results are shown in Table 1. 

Consider at first the upper left-hand corner of Table 1 corresponding to alignment of 
oblate grains when drift tends to be perpendicular to magnetic field 4 (i.e. s < and g < 0). 
For these conditions oj is shown in Fig. 2. The limiting case s — > —0.5, g — > —0.5 gives 

lim (cos 2 ^) = 2 - (17) 
s -> -0.5 

g -> -0.5 

which provides 

5 - 3^2 , , 

a = <jj = ^— « 0.38. (18) 

Note, that for the same values of s and g alignment measure a = 1 was obtained for 
thermally rotating grains in Lazarian (1994a). We would like to emphasis that the 3D plot 
is not symmetrical in respect to interchange of s and g. This reflects the asymmetry in the 
way these two parameters enter the formulae for the alignment measure. 

The lower left corner of Table 1 corresponds to oblate grains tending to drift mostly 
along magnetic field lines (i.e. 5 > and g < 0). The relevant measure of alignment is shown 
in Fig 3. 

It is obvious that s = (i.e. isotropy of external conditions) should produce no 
alignment. We will show this for flakes (g = —0.5). In this case 



(COS 2 0!) 




1 arcsinh y/2s 



2 arctanh 



l + 2s , 



lim^)- 1 ^ - 1/35 3 / 2 - -s)- 1/35 3 / 2 ) = -, (19) 



where series representation of arcsinh x and arctanh x (Granshteyn & Ryzhick 1965, 
1.641 [2] , 1.643[2]) were used. This value of (cos 2 0\) corresponds to the isotropic distribution 



4 When we say "tend to be perpendicular", this means that the direction of the drift is 
within [— arccos l/y/3, arccos l/y/3] if we measure angles from the magnetic field direction. 
Other angles correspond to the drift which in the adopted terminology " tend to be along" 
magnetic field. We adopt this terminology just to give a visual picture corresponding to 
different s. 
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corresponding to the "undamped time" (see section 2). Grain interaction with the flux 
results in J alignment in X\y\Z\ system of reference, which for convenience we will call "gas 
system of reference". In the gas system of reference orientation of J is given by 0\. We 
remind our reader that during "undamped time" J moves in the grain coordinates and 
this motion corresponds to 6 2 changing from to tt. After a crossover J stays fixed in the 
external system of reference but is aligned with the axis of major inertia. For instance, for 
an oblate spheroidal grain this means that the angle between the grain symmetry axis and 
the axis of alignment Zi coincides after the crossover with the angle between J and Z\. As 
the time scale associated with crossovers is considerably shorter that that associated with 
spin-ups, the fraction of grains undergoing crossovers at any given moment is negligible and 
therefore the measure of J alignment for oblate grains 

*J = l(cos 2 9,) - I (15) 

coincides for an ensemble of grain with the Rayleigh reduction factor a (see Greenberg 
1968). Similarly, for prolate grains J gets perpendicular to the grain symmetry axis and 
therefore a = — 0.5<tj. 

The mean value of (cos 2 0\) can be found using W(6i } 62)- Indeed, 
(cos 2 ^) = 

cm rJ 2 rJ 2 ^ 2 (1+ :i^-£)3 /2 - m 

We emphasize, that the alignment happens in two stages, namely, at first J is aligned 
in the gas reference frame (during the crossover), i.e. when there is no alignment of J 
whatsoever in the grain reference frame, then a perfect alignment of J is obtained in the 
grain reference frame (during the spin up). 

Although the alignment is attained over a sequence of short time intervals separated 
by relatively long periods of spin-up this does not make this alignment inefficient. Indeed, 
during crossovers the grain angular momentum is minimal and therefore it is easy to change 
its direction. Moreover, in spite of the fact, that the time scale of crossovers might be 
shorter than the time of grain precession in magnetic field, its lines still represent the axes 
of alignment for grains. Indeed, at any particular moment different grains of an ensemble 
undergo different phases of precession. 

There exist a considerable difference between the way oj is obtained in Lazarian 
(1994a) and the way we find it here. In the former paper it was assumed that the angular 
momentum is directed along the axis of grain major inertia. Evidently this is not true for 
the mechanism discussed above. In fact, crossovers are characterized by a disorientation 
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As the change of grain angular momentum in the course of an individual collision is 
small, it is possible to describe the process of J alignment by the Fokker-Planck equation 
(see Reichl 1980, Roberge et al 1993). As a result, following Dolginov & Mytrophanov 
(1976), we can write 

dj(x,n) 2 dj(x,n) 2 d 2 J(x, n) 

-a^ = £ a, ' (x) ~^ + 4 £ Mx) ~a^' (9) 

where x is a vector in the phase space with coordinates J, cos#i, cos6 2} (fit (fi2 and the 
coefficients a 4 - and presented in the Appendix A. The solution of Eq. (9) is as follows 
(see Dolginov & Mytrophanov 1976) 

/(J, cos #i, cos 6 2} n) = 

consta / J 2 (l+g cos 2 d 2 + s cos 2 ' 



rYn J U+gCOS U2+SCOS »!) / 1n x 

^3/2 «xp ^ 2nb 2 p 2 (l + s+g) J ' 'v 1U / 1 

where s is the external flux anisotropy 

s = -\((p 2 ) - HpI) W) - (pI))- 1 (ii) 

and g is the grain non-sphericity 

</ = ^(« 2 -& 2 )- (12) 

Note that (p 2 ) and (pi) are the averaged squared momentum and its Zi component (see 
Fig. 1) transferred to a grain in an individual collision. Both g and s can vary from — | to 
oo. It is easy to see that g = —0.5 corresponds to flakes and g — > oo to needles. 

Calculations in Dolginov & Mytrophanov (1976) have shown that the solution given 
by Eq. (10) is accurate up to |0.25sg| terms for |s| < 1 and \g\ < 1. The accuracy of the 
solution is of the order s~ 2 for s — > oo when \g\ < and g~ 2 for g — > oo for s < 0. If both g 
and s are large the accuracy is of the order of g~ x or s -1 . To find the distribution function 
for angular momenta W (61,62), one needs to integrate Eq. (10) over the magnitude of 
angular momentum. This integration provides 

W(6 l7 6 2 ) = C(s,g)(l + s cos 2 9 1 +g cos 2 6 2 )~ 3/2 , (13) 

where C(s,g) normalizes the distribution so that 

^, , r^ 2 , n rl 2 , n sin 6 1 sin 6 2 , s 

C(s,g) dOJ A62— 1 2 =1. 14 

Jo Jo (1 + s cos^ 61 + g cos^ 6 2 ) ' 

A peculiarity of "crossover alignment" as compared to the Gold one (see Lazarian 
1994a) is that an elementary process of J orientation happens over a short time interval 
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To summarize, if the velocity of bombarding atoms exceeds (72) 1/,2 ^ff2, which is of the 
order of 2.5 • 10 5 cm s -1 for 72 equal 0.2, the collisions dominate in terms of momentum 
deposited with the grain. In dark clouds, the ratio 72 is much smaller due to the decrease 
of relative abundance of atomic hydrogen and therefore the importance of grain-atomic 
collisions increase. In fact, for sufficiently small 72 ratio, gaseous bombardment rather than 
H 2 formation dominates the random torques. 3 

If inequality given by Eq (8) is not satisfied, it is easy to see that the dependence 
3/2(/c — 1/3) is preserved, but the amplitude of alignment becomes x/(l + x) times smaller, 
where x = (( A J) 2 ex ) / (( A J Zo ) 2 )). If the time between the crossovers is t X} it is possible to 
show that alignment takes place on the time scale of the order of a few t x j arccos(exp(— F)), 
where F is the disorientation parameter introduced in Spitzer & McGlynn (1979) (see also 
Lazarian (1995a) for an explicit expression). 

Accounting for the ISM magnetic field H and grain magnetic moment, which is due to 
the Barnett effect does not alter substantially the picture above. The zero approximation 
would correspond to the angle between H and J being preserved; the changes of it due to 
stochastic torques accounted in the first approximation. 

3. The measure of alignment 

The estimates of the previous section indicate that for rather moderate values of 
gas-grain velocities a corpuscular flux can provide larger increments of angular momentum 
than H 2 formation. Therefore it is natural to assume that H 2 formation causes systematic 
torques, while stochastic torques are mainly due to a corpuscular flux. In other words, our 
model enables us to disregard in the first approximation the stochastic contribution caused 
by H 2 formation during the crossovers. 

The distribution of angular momentum can be characterized by a function f(n } J), 
where n is the number of grain-atomic collisions. In general, the direction of J should be 
defined by angles 61 and (fi in the "gas reference frame" and by 6 2 and <p 2 in the "grain 
reference frame" (see Fig. 1). Henceforth grains will be approximated by spheroids with 
semiaxes a and b. 



3 Moreover, the impact of H 2 formation can decrease considerably if a substantial part of 
H 2 molecules are formed over sites corresponding to aromatic hydrocarbon, which produces 
H 2 molecules with low kinetic energies (see Duley & Williams 1993, Lazarian 1995a for more 
details). 
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where L z is a constant torque along the axis of major inertia. The constant in Eq. (5) is 
equal to the quadratic sum of J x and J y : 

Jl = J 2 X + Jy, (6) 

as this is the value of J Zo at t = (where x and y are the body axes of the grain). Therefore 
the angle between the grain axis and J direction changes from nearly zero to tt during a 
crossover. In other words, during a crossover the direction of J remains the same, while 
the rotation in the grain body axes changes from nearly pure rotation about the axis of 
major inertia to rotation about an orthogonal direction, while the grain gradually flips over. 
Spitzer & McGlynn called this peculiar period "undamped time" as during this time J is 
not clamped to the major axis of inertia. 

Following Spitzer & McGlynn (1979), it is natural to consider dynamical evolution 
given by Eq. (5) as a zero-order solution of the problem, and the dynamical effects caused 
by stochastic character of the applied torques as perturbations of the zero order solution. 
We have seen that the zero-order solution envisages a fixed direction of J in the inertial 
x y z frame, the deviations from which are due to stochastic torques. In the absence of 
external gaseous flux all the directions are equivalent and the deviations are isotropic. We 
will show that in presence of the flux, the deviations become anisotropic and so becomes 
the density of J vectors for an ensemble of grains. 

In the course of a crossover grain angular momentum approaches zero and therefore 
grains become very susceptible to stochastic torques. The first kind of these torques is 
associated with H 2 formation and within the adopted model they are isotropic, i.e. mean 
square increments of the angular momentum ((A J;) 2 ) along axis x 0} y 0} z are equal. On the 
contrary, increments ((A Ji) 2 ) ex ) associated with flux - grain interaction are anisotropic. 

If the flux deposits a fraction k of the squared momentum along the z axis, and one 
may roughly estimate the equilibrium distribution of J vectors using the formulae 



3 ((AJ, ) 2 )+ 72 - 1 ((AJ)L)^ I 
2 3((AJ, )2 )+7 -i ((A J)^) 2 

3 

2 



((AJ)L)(«-l/3) 



(7) 



where 7^ is the ratio of atoms striking the grain to that included in the H 2 molecules 
formed over grain surface. Note, that for pure hydrogen environment 72 = 7i- The 
combination 3/2(/c — 1/3) is the measure of J alignment in the absence of disturbances 
caused by torques induced by H 2 formation, and <r} tends to it if 

<(AJ)L)»7 2 ((AJ, ) 2 ) (8) 
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This means that each grain rotates around its major axis of inertia that we denote z-axis. 
Therefore the problem becomes one dimensional as the two other components of the torque 
contribute only to insignificant nutations. The number of H 2 molecules ejected per second 
from an individual site is ~ ^\l 2 V\riHV~ x , where 71 is the portion of H atoms of number 
density njj approach the grain with velocity v\ is reacted to form H 2 molecules while v is 
the number of active sites over the grain surface. Then the mean square of a residual torque 
is 

([M z ] 2 ) « ^n 2 H m H2 vlEv-\ (1) 

where the coefficient 1 / 4 accounts for the fact that only components of the angular 
momentum parallel to the z axis contribute to M z . The mean squared angular velocity 
yj (0 2 ) of a grain depends on a characteristic time of systematic torques and the frictional 
damping time 

^0.6 C 9 , 2 

o n nmvi 

where S n is the grain cross-section for the gaseous flux,m 3 is the grain mass, m and n are, 
respectively, the mass and the number density of gaseous atoms. In fact, it was shown in 
Purcell (1979, eq.(52)) that 



(02) = ^r-Hutr L ^ (3) 



where I z is the z component of the momentum of inertia. 
If t L > t d : 

= (IM Z ?) 1/2 T- (4) 

Assuming that v ~ 100, 71 ~ 0.2 and E ~ 0.2 eV, one obtains 10 8 s -1 , which 
considerably exceeds the corresponding thermal velocity of rotation. Further on, for brevity, 
we will call "suprathermally rotating grains" just "suprathermal grains". 

2.2. Crossovers 

It was shown in Spitzer & McGlynn (1979), that the direction of angular momentum 
of a suprathermal grain, subjected to a regular torque along its major axis of inertia, does 
not change in axes x y z fixed in space, while its modulus changes as 

J 2 Zo = L 2 f + const, (5) 



- 4 - 



the Gold mechanism in the sense that it is based on depositing of angular momentum by 
a gaseous flux. The difference is that this deposition takes place during short intervals of 
crossovers, when the grain angular momentum is minimal. The efficiency of this process 
depends on the ratio of the kinetic energies of nascent H 2 molecules to the kinetic energy 
of striking atoms and on the accommodation coefficient for atomic hydrogen. In short, if 
stochastic torques are dominated by H 2 formation, the mechanism is suppressed. Contrary 
to this, the second mechanism does not depend on the stochastic torques arising from the 
gaseous bombardment, but on the rate at which atoms arrive at grain surface. We show that 
the frequency of crossovers depend on the orientation of a non-spherical grain in respect 
to a gaseous flux and this causes alignment. It is likely, that for mildly supersonic drags, 
the second mechanism prevail; the opposite is true for hypersonic drags. To summarize, 
our study testifies that suprathermal rotation does not prevent grains from being aligned 
mechanically and we find that high degree of alignment is attainable. 

Below we do not discuss regular torques caused by variations of the accommodation 
coefficient or photoelectric emission. These processes and the alignment that they can cause 
are discussed in Lazarian (1994b). Neither we address the important question of whether 
the mechanical alignment is prevalent within particular regions of the ISM. In our next 
paper in the series we are going to address this problem by comparing the relative efficiency 
of different mechanisms for a number of typical ISM regions. 

The structure of the paper is as follows. In section 2 we briefly discuss the effect of 
gaseous bombardment on a suprathermally rotating grain undergoing a crossover, then in 
section 3 we find analytical expressions for the alignment measure attainable by grain - gas 
interaction during crossovers. Section 4 deals with the alignment caused by grain - gas 
interaction during spin-ups. A joint action of the two mechanisms is discussed in Section 5. 

2. Suprathermal rotation &c crossovers 

2.1. Torques due to H 2 formation 

Suprathermal rotation of the ISM grains was theoretically discovered by Purcell 
(1975). There it was shown that the ejection of H 2 molecules formed over grain active 
sites can cause rotation with suprathermal velocities and this should be the major cause 
of the suprathermal rotation for the ISM grains. Then, in Purcell (1979), it was shown 
that internal dissipation of energy within grains, mainly due to the Barnett relaxation, 
suppresses rotation around any axis but the axis of the greatest inertia on the time-scale 
~ — s, where tj is the ratio of grain rotational energy to the equipartition energy ~ kT. 
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1. Introduction 

Although star light polarization by aligned grains was discovered as far back as 1949 
(see Hilther 1949, Hall 1949), the cause of grain alignment remains still a bit of a mystery 
(see Goodman et. al. 1995). The proposed mechanisms can be subdivided into two classes, 
namely, mechanical and paramagnetic alignment. 1 

At the moment paramagnetic mechanism suggested in its original form by Davis & 
Greenstein (1951) is believed to be more promising. In fact, its suprathermal modification 
proposed by Purcell (1975, 1979) and its superparamagnetic or ferromagnetic modification 
first suggested by Jones & Spitzer (1967) and further developed by Mathis (1986) are 
widely referred to as the likely candidates for explaining the observed large-scale pattern 
of polarization. However, the suprathermal mechanism was shown in Spitzer & McGlynn 
(1979) to produce only a marginal improvement of paramagnetic alignment, while 
superparamagnetic mechanism was critically discussed by Duley (1978). At the same time, 
the original Davis & Greenstein proposal was shown to be inadequate in Jones & Spitzer 
(1967). We believe, that the ambiguous situation with the ISM paramagnetism deserves a 
special discussion which we started in Lazarian (1995a, b) and intend to continue elsewhere. 
The present paper is devoted to mechanical alignment. 

Mechanical alignment pioneered by Gold (1951, 1952) is believed to be not applicable 
to the majority of the ISM grains, which according to Purcell (1975, 1979) should rotate 
with the energies considerably in excess of the thermal one. according to Purcell (1979), the 
dominant reason for this is a quasi-regular torque due to recoils from nascent H 2 molecules 
being formed at catalytic sites over grain surface. Therefore every suprathermally rotating 
grain behaves as a tiny gyroscope which tends to preserve its direction of rotation and 
thus be insensitive to stochastic torques associated with Gold-type processes. 2 However, 
the very fact that the Gold alignment is not applicable does not necessarily entail that no 
mechanical alignment is possible at all. 

In the present paper we consider two processes that can provide mechanical alignment 
of grains, which suprathermal rotating arises from H 2 formation. The first one is similar to 



1 This classification omits a mechanism of ferromagnetic grain alignment suggested in 
Spitzer & Tukey (1951). However it is shown in Lazarian (1994b) that the latter mechanism 
can be only important for dark molecular clouds. 

2 Note that to simplify our presentation we do not speak about Harwit (1970) mechanism, 
which was shown in Purcell & Spitzer (1971) to be inferior to the Gold mechanism for an 
absolute majority of astrophysically interesting situations. 
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ABSTRACT 

It is shown that mechanical alignment of grains can be efficient for grains 
rotating suprathermally, i.e. with kinetic energy substantially exceeding k 
(the Boltzmann constant) over any temperature in the system. The paper 
studies suprathermal rotation caused by H 2 formation and the alignment that 
takes place due to crossover events. Gaseous bombardment in the course of a 
crossover as well as both gaseous friction and poisoning of active sites are shown 
to produce alignment. The first type of alignment happens due to the angular 
momentum deposited by a corpuscular flux with a grain, the second is caused 
by the change of the mean time back to crossover due to the interaction with a 
gaseous flux. We show that the two processes act as to decrease the grain cross 
section in respect to the flux and we find the Rayleigh reduction factor for the 
joint action of the two processes as well as the range of applicability of each 
of the processes. Our study indicates that mechanical alignment can be more 
widely spread than it is generally accepted. 
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